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Abstract. In this work we investigate the evolution of the X-ray emission of a cluster of single young massive 
stars with different metallicities. We have considered the X-ray contribution originated by the diffuse gas heated 
by the mechanical energy released by stellar winds and SN explosions as well as the X-ray contribution from SN 
remnants. The resulting ionizing spectrum (i.e. A < 912 A) has been used to compute the expected intensity of 
the nebular Hell A4686 A. The observational ratio Hell/H/3 could be reproduced by the models assuming that a 
fraction of the mechanical energy produced by the star-formation episode is reprocessed by interaction with the 
ISM as soft X-ray radiation, contributing to the He ionization. However, the discreteness of the stellar populations 
affects the ionizing flux and may be responsible for the observed dispersion of the ratio. We have finally used 
the synthesis models to estimate the contribution of circumnuclear star-forming regions to the multiwavelength 
energy distribution in Active Galactic Nuclei, finding that the UV to soft X-ray continuum in many Seyfert 2 
galaxies seems to be dominated by star-formation processes. 
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1. Introduction 

In the last years convincing observational evidences have 
been collected about the presence of starburst regions in 



or around active Seyfert 2 nuclei (Hcckman et al. 1997; 



Gonzalcz-Dclgado et al. 1998, and references therein). It 
has been found that most of their UV light originates 
in (circum-)nuclear star-formation sites; the possible con- 
nection between the nuclear activity and the properties 
of these starbursts is still a matter of debate. According 
to the unified scheme of Seyfert galaxies, the active nu- 
cleus should be hidden by an opaque torus in the case of 
Seyferts 2, explaining why the collected UV light is domi- 
nated by young, massive stars. The low degree of contami- 
nation by the active source in the UV provides detailed in- 
formation abuul tiro piupeilies uf Llie sLai-fuiinaLiuii piu- 
cesses. Muieuvei, tiAUapulaLing Lu Lhe ladiu - X- ■lay iambus 
the emission associated with the starbursts, it should be 
possible to disentangle the fractional contribution of both 
sources (the starburst and the active nucleus) at differ- 
ent energy ranges. Multiwavelength evolutionary synthe- 
sis models normalized to the observed UV emission would 
be the ideal tool to perform this analysis. 



With these ideas in mind we started some years ago 
a program to extend our evolutionary synthesis models 



( [Arnault et al. 1989| , Paper I, [Mas-Hesse fc Kunth 1991 
Paper II, hereafter MHK, |Cervino fc Mas-Hesse 1994 



Paper III, hereafter CMH) to high-energy ranges (soft and 
hard X-rays, 7-ra ys). The predictions o f 7-ray emission 
were presented in Cervino et al. (2000a| ), and this paper 



is devoted to the X-ray emission. It became evident from 
the beginning that to properly reproduce the X-ray emis- 
sion, both the contribution of evolved binary systems and 
the heating of the diffuse interstellar gas by the release 
of mechanical energy (by stellar winds and supernova ex- 
plosions), had to be considered. Nevertheless, while High 
Mass X-ray binaries contribute essentially to the hard X- 



ray range (above few keV, Van Bever et al. 1999; Van 



Bever & Vanbeveren 2000), the soft X-ray emission is dom- 



chanical energy from the starburst ( 


Heckman et al. 1995 


1996; 


Dickow et al. 1996 




Stevens & Strickland 1998a 


b 


Strickland & Stevens 1999|; Rummers et al. 2001 


). We will 
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therefore discuss in this paper the predicted soft X-ray 
emission, and its effects on other observables, like the rela- 
tive intensity of Hydrogen and Helium emission lines. The 
contribution of binary systems and evolved sources to the 
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hard X-ray emission will be discussed in a forthcoming 
paper. 

Additionally, the statistical dispersion due to the dis- 
creteness of actual stellar populations as presented in 



Cervino et al. (2000b, 2002, 2001) is particularly impor- 



tant in the high energy domain, where the number of effec- 
tive sources is lower. This computed statistical dispersion 
allows to perform a better comparison with real systems 
and to evaluate the statistically relevant sources in each 
energy range. 

We present in Sec. our evolutionary synthesis model 
and show how the X-ray emission and its associated statis- 
tical dispersion have been computed. In Sec. [| we show the 
predictions on soft X-ray emission. In Sec. ^ we explore the 
effects of the soft X-ray contribution on the He n nebular 
emission line. In Sec. |^ we compare our predictions with 
observational data from star-forming and Seyfert galax- 
ies, aiming to disentangle the relative contribution of star- 
bursts to the global energy budget of Seyferts. We finally 
summarize in Sec. ^ our conclusions. 

2. The evolutionary synthesis model and the 
X-ray emission 

We have updated and improved the synthesis model pre- 



sented in 


Arnault ct al. (198£); 


Mas-Hesse & Kunth 


(1991); Cervino & Mas-Hesse (199' 


) with the following 



modifications: 

1. Inclusion of the full set of Geneva evolutionary tracks 



including standard (Bchaller et al. 1992) and enhanced 
mass- loss rates ( |Meynet ct al. 19*94 ). 
Inclusion of metallici ty dependent a tmospher e models 
for normal stars fr om Kurucz (1991 ), CoStar ( [gchaerer 



fc De Kotcr 1997) and the atmo sphere models for WR 
stars from |5chmutz ct al. (1992j ). 



3. Inclusion of an analytical Initial Mass Functio n (IMF) 
formulation using a dynamical mass-bin)^] (see Cervino 



et al. 2001 , and references therein for a full descrip- 
tion). We have also maintained the original Monte 
Carlo formulation. 

Use of parabolic interpolations in time for the compu- 
tation of the isochrones (see |Cervino et al. 2001 , for a 
full study and discussion). 

Evaluation of the dispersion due to the discreteness of 
the stellar populations in all the outputs (see |Cervino 
et al. 2002 and Buzzoni 1989] for more details) 



For this work, we will use only an analytical formu- 
lation of the IMF without Monte Carlo simulations. The 
basic input parameters of the model are: (i) Instantaneous 
burst of star formation, (ii) Salpeter IMF (a = —2.35) 
with a lower mass-limit equal to 2 M Q and upper mass- 
limit equal to 120 M (in) Set of evolutionary tracks with 
standard mass-loss rates from Schaller et al. (1992). 



1 We have used a resolution in the HR diagram of A (log T e g) 
= 0.001 dex and A(logL) = 0.001 dex for alive stars and Am 
not larger than 0.1 and not lower than 2 x 10~ 7 Mq for stars 
that have reached their life-time at the computed age. 



In the following sections we describe how the differ- 
ent contributions to the soft X-ray emission have been 
included in the computations. 

We will consider as soft X-rays all photons with ener- 
gies between 0.07 and 2.4 keV (the ROSAT band), while 
hard X-ray photons will have energies between 2.5 and 10 
keV. We will also use the einstein band (range 0.1 - 3.4 
keV) in Sec. [| and the total X-ray luminosity (0.02-10 
keV), L x . 

2.1. Sources of soft X-ray emission 

The main (persistent) contributors to the X-ray emission 
in a star-forming region with only single-star populations 
will be essentially individual stars, supernova remnants 
(SNR) and hot diffuse gas. In addition, other sources like 
supernova explosions themselves, produce transient peaks 
of X-ray emission. Since these peaks would be of very 
short duration, we have not included them in our calcula- 
tions, but we want to stress that they could provide some 
degree of variability in the high-energy emission of star- 
bursts. Additionally, neither pulsars in young supernova 
remnants nor X-ray binaries have been included in the 
computations, as explained below. 

2.1.1. Individual stars 

In the case of massive stars, strong mass losses produce 
first a cloud of material surrounding the star. In a second 
stage, winds from t he star shocks this material ge nerat- 
ing X-ray emission ( Chlcbowski fc Garmany 1991 ). The 
luminosity in X-rays for stars hotter than Bl may be ap- 
proximated by Chlebowski et al. (1989 ): 



Lx — 10 7 iboi 



(1) 



The maximum contribution from the total population 
of individual stars is about 10 30 erg s _1 Mg 1 at 3 Myr af- 
ter the onset of a burst, assuming an instantaneous burst 
following a Salpeter IMF slope. As we will see, this con- 
tribution is around 2 orders of magnitude lower than the 
emission from other sources and has not been taken into 
account for computing of the total X-ray spectral energy 
distribution. We have instead included the amount of me- 
chanical energy released by the stellar winds in the total 
budget of energy injected into the interstellar medium, as 
discussed later. 

2.1.2. Supernova Remnants 

Supernova explosions (SN) will contribute strongly to the 
whole spectral energy distribution, but only during a very 
short time (the light-curve of a SN decays very quickly in 
few years) . This contribution is furthermore very variable 
in a short time scale and has not been taken into account 
in the models. 

On the other hand, after a SN explosion, there will 
be a remnant of expanding hot gas that will contribute 
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essentially to the X-ray and radio domains during its free- 
fall and adiabatic phases. The time scale of the free-fall 
phase is about 90 yr (Woltjer 1972), and has not been 
taken into account in the models. The adiabatic phase 



has a time scale of (Woltjer 1972) 



2.0 x KTegfj-n "'-yr 



(2) 



where n is the density where the SNR expands and £50 the 
mechanical energy released by the SN explosion in units 
of 10 50 erg. Assuming a mean value of £50 = 10 and n = 1 
atoms cm -3 , a SNR emits in X-rays during 3.4 x 10 4 yr. 
Following MHK, the contribution to the spectral energy 
distribution from SNRs can be approximated by: 



A S n / L v (t)dt = N SN S„ 



(3) 



where JVsn is the SN rate, r the time in which the SNR 
shows X-ray (or radio) emission and S v the emission as- 
sociated to one individual SNR. Eq. | assumes that the 
supernova rate is high enough to maintain a steady emis- 
sion from SNRs during the burst. 

It is important to point out that in the case of the 
analytical-IMF computations, we have used a mass bin 
such that there is, at least, one star suffering a SN explo- 
sion in every step of 10 4 years starting 2 Myr before the 
computed age. This assures that the computation of the 
X-ray emission is correctly made and that there is no bias 
due to the choice of the mass bin. The final output is the 
average value over the used time step. 

For the emission in radio and X-rays of SNR we have 
done the following assumptions concerning S u : 



X-rays: During the adiabatic phase of the SNR evo- 
lution, we have assumed it will emit soft X-rays fol- 
l owing a composite Raym ond-Smith hot-plasma law 
( Raymond fe Smith 1977 ) at different temperatures 
and at a metallicity consistent with the metal abun- 
dance of the region. In accordance with the compila- 
tion of Dickcl et al. (1990 ), we have assumed the same 
average temperature for the hot gas in the cavity of all 
SNR, as given by Eq. [§ 



soft 



0.65/(0.76 keV) v + 
0.175/(0.23 keV) v - 
0.175/(1.29 keV) v 



(4) 



where f(kT) v is a Raymond-Smith spectrum of tem- 
perature kT. We have normalized the emission of each 
individual SNR within the einstein band (0.1-3.4 
keV) to an average value of 10 36 erg s _1 , similar 



to the emissions given for the Cygnus Loop by Cox 
(19Q) and |Chu fc Mac Low (1990Q for SNR in the 
Large Magellanic Cloud. With these values, and as- 
suming that the emissivity within the einstein band 
is roughly constant in the range of temperatures of 
the adiabatic phase, we obtain finally Seinstein = 
3.4 x 10 40 erg s -1 year, which is the value we have 
considered in our calculations. 



We have assumed in addition a harder component with 
kT = 6.4 keV ( |Cox 1972Q originated at the shock front 
of the SNR, accounting for around 20% of the total 
emission. 

— Radio: Following MHK, we have assumed that the ra- 
dio emission originated by a SNR follows a power law 
with index a = -0.9 from 408 MHz to 4885 MHz (6 
cm). We have updated the emission associated with 
a single SNR considered in MHK to a value at 6 cm 
of 3 x 10 9 Jy Kpc 2 year, that takes into account the 
possible contribution from cosmic ray s and interaction 
with other remnants (see MHK and Mas-Hesse 1992 
for more details). 

Individual point sources detected by ROSAT in nearby 
galaxies arc inde ed generally associated with SNR 
( lAschenbach 1998| ). 



2.1.3. Mechanical energy released into the Interstellar 
Medium 



The large amounts of mechanical energy released by the 
starburst process are expected to heat the interstellar 
medium around the newly formed stars. Observations with 
einstein already showed the presence of hot diffuse gas 
within the H 11 shells of 30 Dor ( |Wang fc Helfand 199l| ). 
Combining ROSAT and ASCA data, Wang (1999| ) confirmed 
the thermal nature of the emission, originated by diffuse 
gas at 2 — 9 x 10 6 K, in the soft X-ray band. Similar results 
have been obtained in many other starburst-dominated 
galaxies: Stevens fc Strickland (1998a| p|) found that the 
soft X-ray spectra of a sample of 7 Wolf-Rayet galax- 
ies could typically be well fitted by a thermal Raymond- 
Smith component with temperatures in the range kT = 
0.3-1.0 keV (corresponding to few times 10 6 K). These au- 
thors concluded that a large fraction of the observed X-ray 
emission in these galaxies originated in a hot superbubble 
formed by the combined action of stellar wi nds from the 
massive, young stars in the central cluster. Strickland fc 
[Btcvcns (1999 ) showed with hydrodynamical simulations 
that the observed soft X-ray emission in NGC 5253 was 
consistent with the predictions for young superbubbles 
blown by the starburst. Other starburst galaxies show- 
ing therma l soft X-ray emission are: NGC 1569, with kT 
= 0.8 keV dHeckman et al. 1995| ); Arp 220, kT = 0.61 keV 
(|Hcckman et al. 1996[) or Hcnizc 2-10 with kT = 0.31 keV 



(Dickow et al. 1996) and Mrk 33 = Haro 2 (Summers et 
|al. 2001|) , among others (see also |Heckman et al. 1997| ). 

We have computed the total release of mechanical 
energy from stellar winds following the prescriptions in 



Leitherer et al. (1992). The mechanical energy released by 



SN explosions ha s been taken from the so lar metallicity 
SN simulations of |Woosley fc Weaver (19951) for SN II and 



from t he simulations of Helium bare cores of |Woosley et al 



1995 ) for WR stars. We have subtracted from the energy 



of each SN explosion the total energy emitted in the 0.02 - 
10 keV band by the SNR during the adiabatic phase. It has 
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a value in the einstein band of Seinstein = 1-07 x 10 48 
erg. 

We have also assumed that a variable fraction e* s of 
the total mechanical energy released will end up interact- 
ing with and heating the diffuse interstellar gas to an av- 
erage temperature of 0.5 keV. As a first approximation we 
have thus considered a contribution to the X-ray emission 
from starbursts coming from a Raymond-Smith thermal 
plasma at a fixed temperature kT = 0.5 keV. 

Different estimations of e* s can be found in the litera- 
ture based on hydrodynamical models, but they are quite 



dependent on the assumed input physic (see Strickland & 



Steven ; 199E , 2000, for an extensive review). The standard 
bubble model of Castor et al. (1975 ) , that assumes a con- 
stant energy input into the bubble, shows that only 20% 
of the energy injected into the ISM is used in the bubble 



expansion (i.e. < 0.80). On the other hand Strickland 



& Ste vens (1999[ ) find 6^=0.05 from their hydrodynam- 
ical simulations (assuming also a constant energy input 
into the expanding bubble). 

We have chosen an arbitrary value of e^=0.2 for the 
presentation of this work. But we want to stress that the 
actual value of e* s should be derived by fitting the ob- 
servational values on an object by object basis. It is in- 
teresting to note that a high value of e* s leads to a high 
X-ray emission, but, since the total energy in the burst 
must be conserved, it also leads to a small radius for the 
expanding bubbles. Therefore, the comparison with obser- 
vational data will only be consistent when both the X-ray 
emission and the structure of expanding bubbles are taken 
into account. The effect of changing the efficiency could 
be computed at any time using the relative values given 
in Table |. 

2.1.4. Other X-ray sources 

- High Mass X-ray Binaries (HMXRB): High Mass X- 
ray Binaries become X-ray "active" (permanent) after 
the primary star has collapsed in a compact compan- 
ion, and when the secondary star enters the giant or 
supergiant phases, starting to accrete material onto 
the surface of the compact companion. HMXRB will 
become therefore active only after the first 4-5 Myr of 
evolution, depending on the upper mass limit of the 
IMF. The X-ray emission properties will be dependent 
on the radius of the compact companion. Assuming a 
radius of the compact companion of around 10 km, the 
bulk of the emission will be mainly in the hard X-ray 
band (2.5-10 keV) with a tail that will also emit in soft 
X-rays. The hard X-ray predictions of our models in- 
cluding binary systems and their statistical relevance 
will be presented in a forthcoming paper 

— Pulsars in young SN remnants (YSNR): Van Bever 



fe Vanbeveren (2000) ) show that the loss of rotational 
energy of young pulsars in the form of electromag- 
netic waves are another possible source of X-ray emis- 
sion. This emission can be fitted by a Raymond-Smith 



model with a temperature kT around 2 keV. According 
to these authors, the total X-ray luminosities associ- 
ated to these sources are around 10 32 erg s _1 Mg 1 (for 
a starburst with a mass normalization given by an IMF 
with a slope a = —2.7 and mass range 10 -100 M Q ). 
This emission takes place only after the most massive 
stars have exploded as supernovae, and is negligible 
when compared with the total emission associated to 
both SNR and diffuse hot gas, even for rather low effi- 
ciency values. Furthermore, only a very small number 
of such objects is expected, so that their prediction 
is affected by very strong statistical fluctuations. For 
these reasons we have not considered this contribution 
in this work. 

Low-mass X-ray (LMXRB) binaries and Compton 
scattering: Persic fc Rephacli (2002] ) compute the 
steady-state X-ray emission of star forming galaxies. 
In their work they include the emission produced by 
low-mass X-ray binaries, compton scattering and the 
emission of the compact nucleus in AGNs. They also 
include the contributions from the diffuse emission and 
HMXRB. Since our work aims to establish the contri- 
bution to the soft X-ray from star-forming regions, we 
have not included the emission of compact nuclei. We 
are also restricted to young ages (t < 20 Myr) and sin- 
gle stellar evolution, hence the contribution of LMXRB 
has not been considered. Finally, compton scattering 
is related to the electron density and the presence of 
magnetic fields and will emit mainly in the hard X-ray 



domain (see Persic & Rephaeli 2002, and references 
therein) and has not been included in the current com- 
putations. 



As pointed out before, we have not considered these 
contributions, but we keep in mind that the emission from 
YSNR could be significant in starbursts devoid of inter- 
stellar gas, where no contribution from hot gas is present 
and that HMXRB may increase the emission in both soft 
and hard X-rays. 



2.1.5. X-ray conversion factors 

As far as the X-ray predictions are dependent on the as- 
sumed band, we have summarized in Table [l] the con- 
version factors for the two components discussed in this 
work normalized to the total X-ray luminosity (0.02- 
10 keV). Such conversion factors have been obtained by 
the integration of the assumed X-ray spectrum over the 
given band. The results for all the bands can be found 
in graphic and tabular form in our WWW server at 



tittp : //www. laef f . esa. es/users/mcs/SED, but we will 



show here only the predictions for the X-ray luminosity 
in the ROSAT band, for which most observational data 
are available. In the table we also give the ROSAT bands 
needed for the computation of the ROSAT hardness ratio. 
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Metallicity 






ROSAT 






EINSTEIN 




/ r\ /~i r-r c\ a l i r\ 

(0.07-0.4 keV) 


(0.41-2.4 keV) 


(0.41-1.0 keV) 


(1.0-2.4 keV) 


(0.07-2.4 keV) 


(0.1-3.4 keV) 








SNR 








0.001 


0.290 


0.441 


0.251 


0.190 


0.731 


0.745 


0.004 


0.263 


0.490 


0.301 


0.189 


0.753 


0.759 


0.008 


0.243 


0.526 


0.337 


0.188 


0.769 


0.770 




990 


u.ouo 


0.380 


0.188 


788 


78^ 


0.040 


0.207 


0.591 


0.403 


0.188 


0.798 


0.789 








Diffuse hot gas 






0.001 


0.408 


0.457 


0.360 


0.097 


0.865 


0.813 


0.004 


0.284 


0.606 


0.494 


0.112 


0.891 


0.859 


0.008 


0.221 


0.683 


0.563 


0.120 


0.904 


0.883 


0.020 


0.160 


0.756 


0.629 


0.127 


0.917 


0.905 


0.004 


0.134 


0.788 


0.658 


0.130 


0.922 


0.915 



Table 1. Values of X-ray luminosities in different bands normalized to the total X-ray luminosity (0.02-10 keV) 



2.2. Influence on the ionization structure 

We have also taken into account that X-ray emission 
will produce a relatively small fraction of ionizing pho- 
tons (compared with those generated by massive stars). 
Nevertheless, such photons will be predominantly more 
energetic than the ones of stellar origin, and will affect the 
ionization structure of the nebular gas, and thus the rela- 
tive emission of some lines, as we will discuss below. The 
number of ionizing photons between the ionizing edges of 
H and He are given in Table for the two components 
considered, as obtained from the integration over the as- 
sumed spectrum of each component. In this case, we have 
normalized the values to the total luminosity emitted in 
the ROSAT band. Note that the normalization is different 
to the one used in Table [y. 

As an example for any further use of this table, let us 
assume a Z=0.001 metallicity star forming region where 
the X-ray luminosity produced by SNR is 5 x 10 38 erg s" 1 
and the X-ray luminosity produced by the diffuse hot gas 
is 5 x 10 39 erg s _1 , both in the ROSAT band. There will 
be an additional component in Q(H + ) of 2.07 x 10 48 and 
2.40 x 10 49 photons s _1 due to the X-ray spectrum from 
the SNR and the hot diffuse component respectively. Of 
course, these additional Q(H + ) values are dependent on 
the assumed X-ray spectrum, but it gives us a first order 
approximation of how important the contributions of the 
X-rays sources to the ionizing flux are. 



2.3. Evaluation of the dispersion 

In the following paragraph we summarize how to calcu- 
late the dispersion due to the discreteness of the stellar 
population. We refer to Buzzoni (198S| ) and Cervino et al 



Metallicity 



(2002| , gOOjj ) for further details. 



The IMF gives the probability, w% , of finding a number 
of stars within a given mass range. Each Wi corresponds to 
the mass binning used in analytical-IMF synthesis codes. 
If we assume that each Wi follows a Poissonian distribution 



( Cervino ct al. 2002 ), the variance, of of each Wi is equal 



to the mean value of the distribution, iOj. Let us assume 
now that each star has a property whose contribution to 



0.001 
0.004 
0.008 
0.020 
0.040 



Q(He++) 
(xlO 9 ) 



Q(He+) 
(xlO 9 ) 



Q(H + ) 
(xlO 9 ) 



2.35 
2.24 
2.17 
2.11 
2.05 



SNR 
3.80 
4.05 
4.24 
4.43 
4.53 



4.13 
4.38 
4.36 
4.74 
4.85 



Diffuse hot gas 

0.001 2.65 4.33 4.79 

0.004 2.07 3.50 3.87 

0.008 1.75 3.06 3.41 

0.020 1.49 2.69 2.99 

0.040 1.37 2.51 2.80 

Table 2. Values of the number of ionizing photons at 
different edges normalized to the X-ray luminosity in the 
ROSAT band (0.07 -2.4 keV), as originated by the two 
considered contributions (SNR and diffuse hot gas heated 
by the release of mechanical energy) . 



the integrated property A is given by Widi with a variance 
ufaf — Widf . The total variance of the observable A is the 
sum of all the variances. The relative dispersion is: 



OA 

A 



J2i W * a i 



2\l/2 



1 



VK^A) 



(5) 



where the last term gives us the definition of N e g(A) de- 
scribed by Buzzoni (1989| ) . Note that N e g(A) is normal- 
ized to the total mass. N c g(A) is not a real number of 
stars, but is gives us an idea of how many effective sources 
contribute to any given observable. Cervino ct al. (2002 ) 
show that N e g also defines the mean value of a Poissonian 
distribution that can be used to obtain the correspond- 
ing confidence levels of any observable in function of the 
amount of mass transformed into stars. 

Let us stress that IMF fluctuations are present in 
Nature (the number of stars are always discrete), so, N e g 
is not an evaluation of the errors of the synthesis models. 
It is an evaluation of the dispersion intrinsically present 
in real clusters, i.e. the dispersion is also an observable. 
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This intrinsic dispersion must be taken into account when 
fitting observed quantities to model outputs before estab- 
lishing any conclusion. 

For differential quantities, like the SN rate or the me- 
chanical power, the obtained quantities are the mean value 
averaged over the time step (0.1 Myr in our case). The cor- 
responding dispersion shows the variation over such aver- 
age mean value (Cervino et al. 2001). 



3. Predictions on soft X-ray emission 

As discussed in the previous section, both supernova rem- 
nants and diffuse hot gas will be the main contributors to 
the soft X-ray emission if only single-stars are considered. 
Depending on the evolutionary state of the starburst, one 
or the other component will dominate the emission at a 
given energy. We show in the top panels of Fig. |l]the pre- 
dicted X-ray emission in the ROSAT band as a function of 
metallicity for two cases: (a) X-ray emission produced by 
SNR plus hot gas, with e* s = 20% in the left panel, and 
(b) the X-ray emission produced by SNR only (i.e., with 
e^j = 0%) in the right panel. 

The X-ray emission depends on the star-forming re- 
gion metallicity only during the first few Myr, when the 
emission is dominated by the reprocessing of mechanical 
energy, since stellar winds are strongly metallicity depen- 
dent. On the other hand, after around 5 Myr the effect of 
metallicity on the X-ray emission is small since the super- 
nova rate is essentially independent of metallicity. 

It is clear that the total X-ray intensity is strongly 
dependent on the associated value, as shown in the 
top panel of Fig. ||. It is remarkable that even a relatively 
low value e^ s — 0.05 is enough to produce significant X-ray 
emission, even during the first Myr of a starburst, when no 
other sources are yet active. This is especially interesting 
considering that most star-forming galaxies have been ob- 



serve d at ages between 3 and 7 Myr ( Mas-Hesse fe Kunth 



199£). In general, a change in the reprocessing e^L from 



0.05 to 1 implies approximately one order of magnitude in 
total X-ray emission. 

Since most starbursts seem to have formed massive 
stars according to an IMF with a slope close to Salpeter's 
one ( Mas- Hesse fe Kunth 199^ ), and since their metallic- 
ities can be derived from the analysis of optical emission 
lines, it should be possible to derive a first order estimation 
of the average value by just comparing the predictions 
with the observed soft X-ray luminosities. 

The effective mechanical energy and mechanical power 
remaining available to drive gas flows away are smaller 
than the total amounts generated by the starburst, since 
some fraction of the energy, parameterized by e* s , is re- 
processed into thermal emission, and does not contribute 
to accelerate the gas. If such effect is omitted, the age de- 
duced from kinematical studies will result systematically 
lower than the one obtained from the global analysis of 
the starburst. 

In Fig. |^ we show the mechanical energy and mechan- 
ical power and the corresponding effective mechanical en- 



ergy and power when a value of — 0.2 is used for 
the X-ray emission for a solar metallicity burst. We also 
show these quantities when the correction due to the X-ray 
emission of the SNR is not taken into account (pg. ncorr ) . 

Lower panels of Figs, [l], || and [2] show the value of 
N e g(Lx), N c s(Lk) and N e g(EK). The right axis shows the 
amount of gas transformed into stars for the given IMF 
and mass limits needed to ensure a dispersion lower than 
10%, denoted as M 10% (L X ), M 10% (L K ) and M 10% (E K ). 
We will use this notation for subsequent figures. 

In the case of the X-ray luminosities and the mechani- 
cal power two regimes can be separated. At the beginning 
of the burst the dispersion is dominated by the stellar 
wind component, i.e. there is a large number of effective 
sources that are contributing to the luminosity. When the 
first SN explodes, the value of N e s decreases abruptly: 
as far as the remnant of a single SN produces more X-ray 
emission than the stellar winds associated component, the 
effective number of sources decreases and the dispersion 
due to the discreteness of the stellar population becomes 
larger. It is also consistent with the fact that for evolved 
starburst the possible dispersion will be dominated by the 
occurrence of SN events, and so does the X-ray emission. 
In the case of the mechanical energy, the value of iV c ff (Ek) 
is larger (hence the dispersion lower) because it is an inte- 
grated quantity (see Cervino et al. 2001, for more details). 



3.1. Comparison with other works 

There are at least two other studies predicting the soft 
X-ray emission in starburst galaxies: Bilich et al. (2001] ) 



and |Strickland fe Stevens (1999] , |2000| ). In both cases the 
X-ray emission is obtained by simulations of superbubbles 
and hydrodynamical models taking into account the time 
evolution of the kinetic luminosity. 



In Bilich et al. (2001 ) , the X-ray luminosity is obtained 
from two components: the interior of the superbubble and 
the shell, with an analytical dependence of the metal- 
licity on the radiative cooling function, Ax(Z,T). They 
show that the soft X-ray emission depends on the enrich- 
ment of the hot ISM by stellar winds and SN explosions. 
Ne vertheless, their results are s imilar to the ones obtained 
in Strickland fe Stevens (1999| ). 

The study by |Strickland fe Stevens (1999| ) uses a su- 



perbubble model expanding in a vacuum medium (i.e. an 
external pressure equal to zero) and obtains the X-ray lu- 
minosity from the integration of the density structure and 
the radiative cooling function over the volume of the bub- 
ble. They find a correlation between the mechanical lumi- 
nosity and the soft X-ray emission of Lx(t) ~ 0.05 x Lk(^) 
with a time delay between both quantities (variations in 
L^{t) occurs earlier than variations in Lx(t)). They also 
present an extensive discussion of the factors affecting 
the computed X-ray emissio n, assuming superbubbles (see 
IStrickland fc Stevens (2000| )). Their results (their Fig. 6) 
are quite similar to the ones we have obtained (c.f. Fig. 2 
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Age [Myrs] Age [Myrs] 

Fig. 1. Top panels: X-ray emission in the ROSAT band in erg s _1 Mq 1 as a function of metallicity for an instantaneous 
burst with Salpeter IMF slope in the mass range 2 - 120 M Q for two cases: (a) SNR and equal 20%. (left) and (b) 
only SNR (right). Low panels: iV e ff(Lx) m units of Mq 1 . Right axis show the amount of gas transformed into stars 
needed to obtain a relative dispersion of 10%, M 10 %(Lx)- 



scaled to a 10 6 M burst) with an efficiency of 5%, except 
for the time delay as we have explained above. 

Additionally. 1-D hydrodynamical simulations com- 



puted by Pliischkc (2001) taking into account the envi- 
ronment where the bubble expands (an ambient density 
of 40 atoms cm~ 3 is assumed in their models), show that 
the dissipation of £k(£) is around 80% (which includes not 
only radiative cooling, but also the effects of mass-loading 
inside the bubbles). Note that such effects may change the 
time delay between Lx(t) and L^{t). 

In summary, our approximation parameterized in 
terms of e^ff, although quite simplistic, provides a good 
first order approximation to the more detailed superbub- 
ble simulations. 



3.2. The case of Mrk 33 = Haro 2 

As an example, we have applied our models to Mrk 33. 
It is not the scope of this paper to perform a detailed 
analysis of the properties of this galaxy and it has been 



used only as a first order consistency check. We refer to 
[Bummers et al. (200l[ ) for a more complete study using 
different constraints and synthesis models. 

The first step to apply the model is to obtain an es- 
timation of the age of the burst. For this wc refe r to the 
W(H/3) value used by pVlas-Hesse fc Kunth (1999| ), which 
was obtained through a large aperture and was corrected 
from the contamination by the underlying stellar popula- 
tion. We show in Fig. ^ the predictions of the models as 
a function of age, with the corresponding 90% confidence 
limits for different amounts of gas transformed into stars 
and a metallicity Z=0.008. The predicted W(H/3) values 
have been computed assuming that a fraction of 0.3 of ion- 
izing photons are not absorbed by the gas and the other 
0.7 is transformed in H/3 luminosity following Case B re- 
combination (see Mas- Hesse & Kunth 1999, for details). 

As we can see in the figure, W(H/3) alone constrains 
the age of the starburst to around 4.2 Myr. By using a 
previo us version of our set of models, |m as-Hesse & KunttJ 
(1999|) derived an age around 4.8 Myr. For our case, it 
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Fig. 4. 90% confidence limits of W(H/3) for different 

amounts of gas transformed into stars as a function of 

age, for a cluster with Z=0.008. We have over plottcd the 
3- X-ray emission in the rosat band as a function of observed w(m value of Mrk 33 as given by [M;X> Ur>>( 



, for an myUnlanoouy burst with Salpolei' IMF slope — & Kunth fl99flh 

1 , .1™ nnfnllnnh n, , ,l Hi n i , ,,, i , , 1 1 , , , i AT ^fj ,.\ \ ,, J'' 



and soiar iiielallicily and Us corresponding iVcff(Lx) in 
units of Mg 1 and M 10 %(Lx) values. 



should be enough to constrain the age of the starburst 
within the range 4 to 5 Myr, assuming an instantaneous 
burst. 

The kinetic energy of the neutral expanding gas that 



is pus hed out by the shell in th is galaxy is (see |Lcqucux 
et al. l |99j| |Legrand et al. 1997|) E k (Hl)= 3 x 10 54 r 2 (kpc) 



erg, where r is the size of the expanding shell. Some pos- 



sible values of r have been proposed in the literature: 
(i) Legrand et al. (1997) derived an r value of 1.23 kpc 



from the expanding H a shell an d using a distance of 19.5 
Mpc (ii) Bummers et al. (2001 ) found a extended X-ray 
source with an ellipsoidal shape and diameters of of 2.3 
kpc and 1.9 kpc assuming a dis tance of 22 Mpc (i.e . 
r = \fab — l.Okpc) and (iii). Also Bummers et al. (2001 ) 
obtain an average radius 1.1 kpc from B band images. 
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From this last study, we assume the same distance and a 
value of 1.1 kpc and so a E k (Hl)= 3.63 x 10 54 erg. 

Finally, the observed soft X-ray luminosity in the 
ROSAT band ranges from 0.2 to 1.4 xlO 40 erg s _1 . The first 



value corresponds to the HRI observations of Summers ct| 
al. (20()l]), while the second one was derived by Strickland 



& Stev |ens (1999| ) from PSPC data. 

Our model predicts Ek between 3.4 x 10 48 and 5.4 x 
10 48 erg Mq 1 at 4 and 5 Myr respectively, and Lk between 



7.0 x 10 34 (at 4 Myr) and 4.5 x 10 34 erg s^ 1 M" 1 (at 5 
Myr). The uncertainty in Lk is lower than 5% in a 90% 
confidence level and lower than 2% (also in the 90% con- 
fidence level) for Ek assuming a mass transformed into 
stars, Mt ra ns, larger than 10 6 M using a Salpeter IMF 
slope with mass limits from 2 to 120 M 



lines have been proposed as age indicators, like [O in] ones 
by Stasinska & Lcithercr (199E). Additionally, the strong 
soft X-ray flux produced by the hot diffuse gas and SNR 
will contribute to harden the overall ionizing continuum, 
potentially affecting the ratios of different emission lines. 

We show in Fig. [sjthe evolution with time of W(H/3) as 
a function of the metallicity for a model where X-ray emis- 
sion (neither SNR nor mechanical energy reprocessed in 
X-rays) is not taken into account (standard) and a model 
with a value of e* s = 0.2. It can be seen that the presence 
of hot gas may affect the predicted W(H/3) values after 
the first 6 Myr of evolution, i.e., when the ionizing contin- 
uum of the massive stars decreases^]. A first implication of 
these results is that W(H/3) becomes a very uncertain age 
indicator for values below about 10 A. A similar conclu- 



A first order estimation of and the mass of the sion has been obtained by Van Bever et al. (1999) based 



cluster can be obtained using the following relations: 



e eff 



Eobs 
K 



E 



model 



M t , 



T obs 



f mode 



x M tr8 



(G) 



Note that, for a fixed value of M trans , larger values of 
Ek produce larger values in e^ s , and lower values in Lk 
produce larger values in e* s . 

Assuming a constant e^j value along the evolution, it 
is possible to obtain e* s values between 0.02 (at 4 Myr 
with HRI data and Mtrans = 1-1 x 10 6 M ) and 0.31 (at 
5 Myr with PSPC data and M trans = 0.99 x 10 6 M s ). On 
the other hand, M tran s ranges from 0.7 x 10 6 M (at 5 
Myr with HRI data and = 0.04) to 1.3 x 10 6 M (at 

4 Myr with PSPC data and e* g = 0.15). 

So, the observed values can be explained with a 4- 

5 Myr old burst with a mass transformed into stars of 
0.7-1.2 x 10 6 M (following a Salpeter IMF slope in the 
mass range 2 - 120 M ) and a value from 0.02 to 0.31. 

In this simple approximation we have not considered 
observational errors on Ek, neither different treatments of 
the mass-loss rates in the evolutionary tracks nor differ- 
ent velocity laws in the computation of Ek and Lk, which 
would presumably increase the range of possible values. 
So, a more careful analysis, based on a self-consistent com- 
parison of the multiwavelength spectral energy distribution 
and a complete set of observables of Mrk 33 with the model 
predictions is necessary to derive any firmer conclusion. 



4. Predictions on optical emission lines intensities: 
H/3 and Hen 

The H/3 equivalent width has been used by many authors 
in the last years as a reliable indicator of the evolution- 
ary status of a starburst, since it relates the most massive 
stellar population (main contributors to the gas ioniza- 
tion) with the total population in the cluster, which p ro- 
duces the optical continu um (see |Copetti et al. 19*86 or 
Mas- Hesse fc Kunth 1999 as an example). Other emission 



only on the effect of binary systems, but in their case, the 
additional ionizing flux comes both from the donor star of 
the binary system, that becomes (depending on the mass- 
transfer scenario) a WR star, and from the gainer star, 
that becomes more massive hence hotter than before the 
mass-transfer episode. They obtain a higher value for a 
reliable use of W(H/3) (i.e. not dominated by binary sys- 
tems) around 75 A. 

Additionally, the dispersion of W(H/3) when the X- 
ray emission is taken into account becomes larger than 
the one without such hot gas emission. It is due to the 
intrinsic uncertainty of the X-ray emission itself (i.e. the 
small number of effective sources). It also means that the 
age determination of systems with low W(H/3) values is an 
intrinsically difficult task, as far as the possible ionizing 
sources are not only massive stars, but also diffuse hot gas. 

While the inclusion of the soft X-ray flux in the ion- 
izing budget does not affect only weakly the total inten- 
sity of the H/3 line, it can affect more significantly other 
emission lines with higher ionization potentials, like [O in] 
A5007 A. At 54.9 eV, the soft X-ray flux might become 
the dominant contribution to the ionizing continuum, af- 
fecting significantly the expected intensity of the line. 

The nebular Hen 4686 A emission line has been de- 
tected in few Star Forming Galaxies (SFG) ( ponti 1991 



[Schacrcr et al. 1999 ), but is apparently absent in most 
of them. The average ratio of the observed intensity of 
Hen over the inte nsity of H/3, I(He n)/I(H/3), is around 
0.02 ( French 1980 ). Massive stars are not hot enough to 
produce the hard ionizing continuum required for such rel- 
atively large rati os, except perhaps during the WR phase, 
as discussed by Schaerer fc Vacca (1998| ). Nevertheless, 



no clear correlation has been found between the detection 
of the line and the presence of WR stars in the region. 
Additional, harder contributions to the ionizing contin- 
uum are therefore required to explain the observed ratios. 



2 The evaluation of the differences in W(H/3) between models 
that include X-ray emission and standard ones depends on the 
metallicit y and e%r values. We refer to the data in ou r WWW 
server at http: //www . laef f . esa. es/users/mcs/SEE for pre- 
cise comparisons. 




Age [Myrs] Age [Myrs] 

Fig. 5. W(H/3) evolution as a function of the metallicity for a model where the X-ray emission is not taken into account 
and a model with a value of = 0.2 and its corresponding iV e ff(W(H/3)) in units of Mg 1 and M 10 %(W(H/3)) values. 



We have explored the effects associated with the repro- 
cessing of mechanical energy in the interstellar medium 
and the X-ray emission of SNR. The results are shown in 
Fig. |. 

The ionization by hard photons produced by the dif- 
fuse gas heated by the release of mechanical energy into 
the medium can lead to significantly higher values of the 
I(Hen)/I(H/3) ratio sp ecially for low metallicit y clusters. 
As has been shown by Schaerer fc Vacca (1998 ), for solar 
metallicity starburst, a value of the ratio of 0.01 can be 
explained only with the presence of WR stars, but such 
value is never reached in low metallicity clusters. We show 
here that it is possible to obtain such ratio if the soft X-ray 
emission is included. 

However there are some points we would like to note 
here: 

The extension of the ionized region emitting the He n 
4686 A line is generally much smaller than the region 
ove r which H /3 is emitted, as shown by Mai'z-Apcllaniz 
in NGC 4214. To properly compare with 



1. (1 



4686 A flux has to be divided by the H/3 intensity in- 
tegrated over the whole ionized area. Such integrated 
ratios can be smaller by almost an order of magnitude 
than the ratios measured just locally. Another illustra- 



tive example can be found in Luridiana et al. (1999) 



et 

the predictions of synthesis models, the observed He n 



where the authors apply a photoionization model to 
NGC 2363 with the aim of explain with a single model 
all the observations of the region obtained by other au- 
thors using different apertures. Whereas the observed 
I(Heii)/I(H/3) ratios varies from 0.010 to 0.004, the 
global I(Hen)/I(H/?) ratio of the model that fit all the 
observations has a value of 0.001. 
Our models compute the predicted soft X-ray emission 
assuming the released mechanical energy is reprocessed 
immediately into high energy photons. A more realis- 
tic scenario would imply that the mechanical energy is 
first accumulated in the form of energetic gas outflows, 
being released in the form of soft X-rays only when this 
expanding gas interacts with static interstellar clouds. 
Therefore, at some stages, the soft X-ray contributions 
to the ionizing continuum could be stronger than val- 
ues predicted by our models. 
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Fig. 6. I(He n) /I(H/3) ratio as a function of age and metallicity for a model where the X-ray emission is not taken 
into account and a model with a value of e^g = 0.2 and its corresponding Af e ff(I(Hell)/I(H/3)) in units of Mg 1 and 
Mi %(I(Heii)/I(H/3)) values. 



The hard ionizing flux intrinsically has a big disper- 
sion. Its value is controlled by a small number of ef- 
fective sources (WR stars, hot gas clouds or binary 
systems), hence suffers from severe statistic problems. 



In order to support the last items, we show in 
Fig. ^ the 90% confidence limits obtained from the 
AT ff(I(Hen)/I(H/3)) value for a cluster with a metallic- 
ity Z— 0.001 and different amount of gas transformed into 
stars. 

We therefore conclude that the additional hard ioniz- 
ing photons produced by the diffuse gas can explain, at 
least partially, the intensities of Hen 4686 A observed 
in low metallicity star forming regions. However, at this 
point it is impossible to establish the source(s) of the He n 
4686 A. From a theoretical point of view, it is necessary 
to include other sources of X-ray emission that were not 
considered in this work, and establish their statistical rel- 
evance. From the observational point of view, more ob- 
servations in different wavelengths domains are needed to 
constrain the sampling effects due to the discreteness of 



the stellar population. They must be unaffected by addi- 
tional aperture effect biases. 



5. Star- formation processes in Seyfert galaxies 

We have compared the predictions of our models with the 
multiwavelength data for a sample of AG Ns (QSO, Seyfert 



1 and Seyfert 2) and SFGs compiled by Mas-Hesse et al 
|(1995| ). We show in Fig. || the predicted evolution of the 
Lv(0.1 - 3.5 keV) over z/L„(1450 A) ratio for £^=0.2 and 
different mctallicities. The model predictions correspond 
to a Salpeter IMF. We have also plotted the highest (0.040 
metallicity with e* s =l) and the lowest (0.001 metallicity 
with e^f=0 including the SNR component only) values of 
the ratio, and the corresponding 90%CL for a 10 5 M Q clus- 
ter, so that, for a given age, both solid lines (and bands) 
represent the upper and lower limits we should expect as- 
sociated to_a^tar-fc«nadon_episode. Data points are from 
SFG in Mas-Hesse et al. (1995| ) compilation. We have de- 
termined the ages of the plotted objects from the W(H/3) 
values taken from the literature, but it can be seen that 
in any case the results are weakly dependent on the as- 
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90% CL, Z=0.001 

10 4 Mo 




Age [Myr] 

Fig. 7. 90% confidence levels of the I(Hen)/I(H/3) ratio 
as a function of age for a Z=0.001 model with a value of 
= 0.2 for a cluster where 10 4 , 10 5 and 10 6 M© have 
been transformed into stars. We also show the mean value. 

sumed age. UV data have not been corrected for internal 
extinction, so that the plotted ratios have to be taken 
as upper values and should be indeed intri nsically lower 



Age (Myrs) 

Fig. 8. Ratio L„(0.1 - 3.5 keV) over i/L„(1450 A) com- 
puted for an efficiency in the reprocessing of mechanical 
energy into soft X-rays emission of e* s =0.2 and different 
metallicities (with lines as in Fig. ^). In grey, the 90% CL 
for a cluster where 10 5 M Q have been transformed into 
stars for two extreme metallicities and e*„ values. Data 



by fac tors be tween 2 and 5, approximately ( [Mas-Hesse fc 
Kunth |l999|) 



From Fig. || we see that almost all galaxies fall within 
the predictions of synthesis models. It is important to 
stress that some conversion of mechanical energy into 
X-rays is always required. Otherwise the models would 
severely underestimate the observed X-ray/UV ratios. 
Additionally, some galaxies fall outside the model limits. 
It can be due to statistical effects and/or the presence 
of binary systems that will increase the X-ray emission 
and/or the underestimation of the observed UV flux due 
to extinction effects. We want to point out that the effect 
of the extinction in our sample has not been considered. 
Even a moderate extinction would significantly affect the 
UV continuum, putting so the X-ray/UV ratio within the 
model limits without additional X-ray sources. 

Note that our models compute only the persistent 
X-ray emission form single-stellar populations. Transient 
sources like Be/X-ray binaries and SN explosions may in- 
crease the X-ray emission in these galaxies, but only dur- 
ing very short periods of time. In the case of Be/X-ray 
binaries the component may be highly variable and would 
affect only the soft X-ray emission. 

In order to disentangle the relative contribution of 
(circum-)nuclear star-formation processes to the total en- 
ergy budget of low activity AGNs, we have compared their 
L l/ (0.1 - 3.5 keV) over ^L l/ (1450 A) ratios with the pre- 
dictions of our synthesis models. We show in Fig. || the 



points are from SFG in Mas-Hesse et al. (1995) compila- 
tion. Note that the UV continuum of the data points has 
not been corrected from internal extinction. 



corresponding histograms for the Sey 1, Sey 2 and SFGs 
taken from the Mas-Hesse ct al. (1995| ) sample. The ver- 
tical line corresponds to the maximum value predicted by 
the models, as shown in Fig. ||. 

The segregation of Sey 1 and star-forming galaxies be- 
comes evident at first glance: while most of the SFGs in the 
sample fall below the predicted upper value, all Seyfert 1 
galaxies show clearly higher values, by one to two orders 
of magnitude. Therefore, while the physical processes as- 
sumed to take place in star-forming galaxies may explain 
their observational properties, other high energy sources 
are required in the case of Seyferts 1, as expected. 

The case for Seyfert 2 galaxies is more interesting: 
while some of them are located on the plot in the region 
associated to starbursts, some others show a clear X-ray 
excess with respect to star-formation dominated objects. 
The presence of star-formation regions in or around the 
nuclei of some Seyfert 2 galaxies has become evident in 
the last years. Hcckman ct al. (1997] ) showed that around 
70% of the UV continuum in some Seyfert 2 galaxies were 
contributed by young, massive stars, with only a small 
contribution being associated to the nuclear active source. 
|Colina ct al. (1997|) estimated that only 1-10% of the UV 



flux in 4 Seyfert 2 galaxies imaged with the HST was orig- 
inated by the nuclear source, while the rest was due to 
young hot stars recently formed around the nuclei. The 
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log [Z„(0. 1-3.5 keV)/^ J/ (1450 K)] 

Fig. 9. Histograms for the ratio L y (0.1 - 3.5 keV) over 
^L^(1450 A) for different types of emission line galaxies. 
Data taken from Mas-Hesse et al. (1995). The UV con- 
tinuum of the galaxy sample has not been corrected from 
internal extinction. The solid line corresponds to the up- 
per value predicted by our evolutionary synthesis models, 
as explained in the text. 



contribution of star-forming regions to the budget energy 
dis tribution of AGNs was also pointed out some years ago 
by Mas-Hesse et al. (1995 ) based in mul tivawelcngth anal 
ysis of these ty pe of galaxies, and by Cid-Fernandes fc 
Terlevi |ch (19951 ). 

Our results show that not only the UV continuum 
seems to be strongly contaminated by the emission as- 
sociated to starburst episodes in Seyfert 2 galaxies, but 
also their soft X-ray emission. This result is consistent 
with the scenario assumed for Seyfert 2 galaxies, since 
the opaque torus surrounding the nucleus would hide not 
only the UV emission associated to the active source but 
also a significant part of its soft X-ray emission due to 
absorption by neutral Hydrogen. For example, a detailed 
analysis of the UV and X-ray emission of NGC 1068 shows 
that a significant fraction of its soft X-ray emission could 
originate from the starburst episodes present around its 
nucleus. On the other hand, the hard X-ray flux would be 
underestimated by our models by more than two orders 
of magnitude, indicating that it is mostly contributed by 
the active source (Jimenez-Baifon et al., in preparation). 
A more extended analysis of some Seyfert 2 galaxies will 
be presented elsewhere. 



6. Conclusions 

In this work we have explored the X-ray emission orig- 
inated in a star-forming region with only single stellar 
populations. The mechanical energy injected into the in- 
terstellar medium by stellar winds and supernova explo- 
sions will heat the diffuse gas to very high temperatures, 
and will be finally reprocessed into soft X-ray emission. We 
have found that the X-ray emission observed in starburst 
galaxies can be well explained assuming that a moderate 
fraction of the mechanical energy is finally reprocessed 
into X-ray emission. The rest of the mechanical energy re- 
leased leads to the expansion of the gas, creating bubbles 
and gas flows at galactic scales. An interesting implication 
is that the age deduced from kinematical studies will re- 
sult systematically lower than the one obtained from the 
analysis of the emission line spectrum if the heating of 
the gas and the X-ray emission is not taken into account. 
Additionally, the inclusion of X-ray transient systems, like 
SN explosions and Be/X-ray binaries showing bursts of X- 
ray emission, may lead to some degree of variability in the 
high energy emission of starburst galaxies. 

We have computed the expected intensity of the neb- 
ular Hen A4686 A, including the additional ionization of 
the gas by the soft X-ray emission originated in the diffuse 
gas. The observational values of the He n/H/3 ratios can be 
reproduced by our models assuming moderate efficiencies 
(about 20%) in the reprocessing of mechanical energies 
into X-ray emission. Alternatively, the mechanical energy 
reprocessing could not be a continuous process. This en- 
ergy released by the massive stars could be "accumulated" 
in the form of accelerated gas flows, which would release 
all this accumulated energy only when the outflowing gas 
interacts with the static interstellar medium. Under this 
scenario the soft X-ray emission originated by the shocked 
gas would provide enough additional ionizing power to ex- 
plain the observed Heii/H/3 ratios. We have also shown 
that such ratios have a high intrinsic dispersion and that 
a deeper statistical study is necessary to investigate the 
source of the nebular He n A4686 A line. 

We have compared the predicted soft X-ray emission 
with observational values for a sample of star-forming and 
Seyfert galaxies, aiming to disentangle the contribution of 
star-formation episodes to the total energy budget of low 
activity galaxies. We have found that while the high en- 
ergy emission of Seyfert 1 galaxies is clearly above the pre- 
dictions of starburst models, both the UV and soft X-ray 
emissions of many Seyfert 2 galaxies are apparently asso- 
ciated mostly to the (circum-)nuclear starburst episodes 
known to be present in these objects. The active source in 
Seyfert 2 galaxies would dominate therefore only at harder 
X-rays. 
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